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In this paper, we address the problem of designing a quantum encoder that maximizes the min- 
imum output purity of a given decohering channel, where the minimum is taken over all possible 
pure inputs. This problem is cast as a max-min optimization problem with a rank constraint on 
an appropriately defined matrix variable. The problem is computationally very hard because it is 
non-convex with respect to both the objective function (output purity) and the rank constraint. 
Despite this difficulty, we provide a tractable computational algorithm that produces the exact op- 
timal solution for codespace of dimension two. Moreover, this algorithm is easily extended to cover 
the general class of codespaces, in which case the solution is suboptimal in the sense that the subop- 
timized output purity serves as a lower bound of the exact optimal purity. The algorithm consists of 
a sequence of semidefinite programmings and can be performed easily. Two typical quantum error 
channels are investigated to illustrate the effectiveness of our method. 

PACS numbers: 03.67.Pp, 02.60.Pn 



I. INTRODUCTION 

The efficient transmission of quantum states over a 
noisy channel is a central subject in quantum informa- 
tion technologies The mathematical description of a 
quantum input-output relation is as follows. Let Ti and /C 
be finite-dimensional Hilbert spaces of an input quantum 
state and the corresponding output, respectively. We de- 
note by C{Ti, K.) the set of linear operators from H to K,, 
and S{'H) the set of quantum states on H. The Marko- 
vian evolution of a quantum state p G S{Ti.) through a 
quantum channel A is typically modeled using the Kraus 
representation as 



P 



(1) 



where the Kraus operators Ai e C{H,JC) satisfy 
J2i = In with I-H denoting the identity operator on 
Ti. The purity of a state p is defined as p[p] := Tr (p^), 
which is equal to one if and only if p is pure. Due to the 
decoherence caused by A, a pure input state p = \(f>){4>\ 
may be transmitted to a non-pure output p' = A{\4'){(j)\) 
with p[p'] < 1. It is considered that p{p'] quantifies an 
intrinsic measure of the amount of decoherence induced 
by the error channel A. In particular, this paper focuses 
on the optimal purity: 

(2) 



P{A) max min Tr [A{m{^,\f 



where the minimization with respect to the state |0c) 
takes into account the worst-case scenario of informa- 
tion processing. The maximization with respect to the 
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codespace C G Ti. is motivated by the fact that we often 
have an opportunity to decrease the effect of decoherence 
by encoding our information into a higher-dimensional 
space; this is suggested by the theories of quantum error 
correction (QEC) [l|, i, i, i] and decoherence-free sub- 
space (DFS) [1, Sill- For example, embedding an input 
state = 01 10) + (j)2\l) & into a codespace spanned 
by 1 00) and |11) through the encoding process 



,100) 



|11)gCc7^ = C'' (3) 



appears to improve the output purity. Clearly, the most 
desirable situation is the existence of a DFS, i.e., a 
codespace that satisfies P{A) — 1; but unfortunately this 
is a rare case. In this sense, the optimal codespace C is 
regarded as the best possible approximation of a DFS. 

However, the max-min problem ^ is very hard to 
solve because it is non-convex with respect to both C 
and \(j)c)- To understand the structure of P{A), in Q 
Zanardi and Lidar considered channel purity for a fixed 
codespace C as 



PiA,C) 



min Tr [^(|0e) (0c|)' 

\<Pa)ec 



(4) 
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and derived the alternative expression 

P{A,C) = min (0,1 (g) {<j>,\n{A)\<l>c) ® |0c), 

where the Hermitian operator Cl(A) is defined by 

n{A) := Y.(A]A,) ® (AjAj) G C{U^^,U®^). (5) 

This expression was used to derive a bound on P{A,C) 
in terms of U,{A) and C, using techniques to calculate the 
expectation value of the "Hamiltonian" Q.{A). In the spe- 
cial case where eigenvectors of Vl{A) are product states in 



2 



a symmetric subspace of analytical expressions for 
P{A,C) were obtained. However, in general the max-min 
problem ([2]) does not have an analytical solution, leading 
us to take a computational approach. 

From a computational point of view, owing to the rapid 
progress of computers, there have been many recent ad- 
vances with a great potential for solving important prob- 
lems in quantum theory. Convex optimization, and in 
particular semidefinite programming (SDP) [lalUl.have 
proven useful for quantum optimization problems such 
as a test for distinguishing an entangled from a separa- 
ble quantum state [H, Uj, [13, [IE and a design of 
optimal measurement in linear quantum systems p,7i | . In 
addition, in [l^, [l^, [S^l some quantum error-correction 
problems were solved using SDP, taking advantage of the 
well-known convexity of a set of quantum channels known 
as the Jamiolkowski isomorphism pll |. 

In this paper, we first use the same convexity property 
to set up a non-convex optimization problem that cap- 
tures our goal and all the constraints. Then, we provide 
an algorithm that computes an exact local optimal solu- 
tion of the hard non-convex problem ^ for the codespace 
of dimC = 2. This implies that the exact global optimal 
solution of ^ can be obtained by appropriately choos- 
ing an initial condition of the algorithm. The algorithm 
is represented by an iterative SDP and is thus computa- 
tionally tractable. The derivation of the SDP consists of 
two stages. The first one transforms the constraints to 
equivalent Linear Matrix Inequality (LMI) constraints. 
The key idea used to obtain the LMI in this stage is 
the Sum- of- Squares characterization of a polynomial con- 
straint [13, [23, HI]- In the second stage, a non-convex 
rank constraint of the matrix variable is tackled via the 
log-det (logarithm of determinant) heuristic psi. I26ll27t. 
Furthermore, we will show an extended version of the 
above SDP algorithm that computes a lower bound of 
the optimal purity P{A) for the general class of C. 

This paper is organized as follows. Section II reviews 
the Jamiolkowski isomorphism, which is used to formu- 
late the optimization problem in Section III. The SDP 
algorithm is presented in Section IV. The general case 
that leads to a suboptimal solution is discussed in Sec- 
tion V. In Section VI, we examine two typical quantum 
error channels, the hit-flip channel and the amplitude 
damping channel, and demonstrate the effectiveness of 
our method. Section VII concludes the paper. 

Notation: A Hermitian matrix X = e /:(C",C") 
is positive semidefinite if (a|X|a) > 0, V|a) G C"; the 
inequality X >Q represents the positive semidefiniteness 
of X. We use /„ to denote the n x n identity matrix, 
which is the same as lu when dimTi = n. For a matrix 
X = (xij), the symbols X^ and X* represent the ma- 
trix transpose and the elementwise complex conjugate of 
X, i.e., X^ = (xj,) and X* = {x*j) = {X^^ , respec- 
tively; these rules are applied to any rectangular matrix 
including column and row vectors. 5R(X) and ^{X) de- 
note the real and imaginary part of X, respectively, i.e.. 



II. THE JAMIOLKOWSKI ISOMORPHISM 

The main purpose of this section is to review the fol- 
lowing important fact known as the Jamiolkowski iso- 
morphism [2l| ; the set of all finite-dimensional quantum 
channels has a one-to-one correspondence with a convex 
set of positive semidefinite matrices ac ting on /C (8) 7i. 
This fact can be seen in various ways [3^ [H, [2^, [sO] ■ 
Here we follow the notations in [3, liSl and obtain two 
matrix representations of a quantum channel, which we 
later use to set up the optimization problem. At the end 
of this section, we present a characterization of quantum 
channels that preserve pure states. 

We consider a general trace-preserving quantum chan- 
nel that maps an input p E S{H) = 5(C") to the output 



(6) 



Let {| j)}i=i,- - ,n and ,m be orthonormal bases 

in H and /C, respectively. Then, any vectors in 7i®^ 



and /C*^ are expressed as |<I>)) = j=i 



\j) and 



1$')) = J2Tj=i'f'ij\^) ® b)' respectively. We sometimes 
use \i)\j) as a short-hand for \i) (g) Let us now define 
the following two specific vectors: 



rn 



(7) 



(8) 



These vectors have the property of being independent of 
the selection of orthonormal basis; for any two orthonor- 
mal bases {\ai)} and {\hi)} in H, we have 



\In)) = Yl ® = 1^^) ® 



(9) 



Note that the invariant property © is not satisfied if 
\Iu)) is defined without the complex conjugation. The 
vectors ([7]) and are related by 

{X®In)\In)) = {lK®X^)\lK)), yXeC{n,IC). (10) 

Further, the following equation holds: 

{{In\{X®In)\lH)) =Ti-X, ^iX e C{n,n). (11) 

We now define a positive semidefinite matrix Xi associ- 
ated with the Kraus operators Xi G C{'H,JC) as 

Xi := E(X, (g> lH)\In)){{In\{X^ (g> Ih)^ 

i 

Then, the trace-preserving condition ^ - xlXi — I-h cor- 
responds to Tr K.Xi = I-u, and the quantum channel ^ 
is expressed in terms oi Xi as 



mX)h = (^u +<-)/2 and (3(X))„- = (x. 



,)/2i. 



' = Trw[(/K;0/)Xi 



(12) 
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Conversely, it is known that any positive semidefinite ma- 
trix Xi g £(/C ®T-L,JC®TL) corresponds to a quantum 
channel with input-output relation given in Eq. (fT2|) . 
That is, there exists a one-to-one correspondence between 
a quantum channel from 7i to /C and a positive semidef- 
inite matrix on /C 7i. 

We next introduce another matrix representation of 
the quantum channel, which will be denoted by X2- To 
this end, we define a vector associated with a quantum 
state p e S{T-L) as 

\p))--^{p®In)\In))^'H^^- (13) 

The vector |p)) is obviously in one-to-one correspondence 
with p. In particular, from Eq. ([9]), the vector represen- 
tation of a pure state p = \a){a\ is given by 

\p)) = {\a){a\ ® Ih)Y.\') ® N)* = l«) ® !«)*• (14) 

i 

In addition, the purity p[p] — Tr (p^) is simply the 
squared Euclidean norm of \p)): 

p[p]=Trip^) = {{p\p)), (15) 

due to Eq. pip . Thus, a quantum state \p)) is pure if and 
only if {{p\p)) = 1. Let us now define X2. Multiplying 
\I)c)) on both sides of Eq. ©, we have {p' ® Ik)\Ik)) — 
Yjii^iP^i ® Ik)\Ik)), which is rewritten by 

[p' ® Ik)\Ik)) = Y.^X, (g> I^)ip (g> I^){In ® X*)\In)) 

i 

^Y.{X,®X*){p®In)\In)), 

i 

because of the property ([T(7)l . Hence, defining the matrix 

X2 ■.= Y,x^<E)X* e /:(H®^/c®2), 

i 

the quantum channel ^ is represented by 

H'^'3\p))^\p'))=X2\p))eIC^'. 
The trace-preserving condition is then given by 

((I,C\X2 - 51 ((^^K^' ® ^^)(^« ^**) 

i 

= Y.^{In\{In ® Xj){In ® X*) = ((/«|. 

i 

The matrix X2 is related to Xi through the following 
rearrangement rule of the matrix elements: 

m\*x2\mY = m\*xsm*- 

This relation is independent of the selection of and 
{|i)} due to Eq. As the rearrangement map is ob- 

viously linear and homeomorphic, Xi and X2 have a 
one-to-one correspondence with each other. We denote 



this relation by Xi = $(^2). The above discussion is 
summarized as follows. 

Lemma 1. Any finite-dimensional quantum channel 
from H to /C is represented by H®'^ 3 \p)) \p')) = 
X\p)) <^ K.^^, where X is in the set 

X{n,lC) = {x e £(H®2,/C®2) 

$(X)>0, ({I^\X = {{In\ }. 

The linear transformation ^{X) is defined with respect 
to orthonormal bases {\i)} € H and {|i)} € /C as 

m*x\k)\er = mrwrm*- 

Clearly, X{Ti,JC) is a convex set with dimension 
m?ii? — n?. It should be noted that a cascade connection 
of two quantum channels X g XiTi^K.) and Y g A'(/C, V) 
is simply represented by the multiplication of those ma- 
trices: Yx e x{n,v). 

Finally, we provide a characterization of quantum 
channels that preserve pure states, i.e., p[p] — p[p'] — 1, 
as follows. 

Lemma 2. For a quantum channel X E X{H,IC), the 
following three conditions are equivalent. 

(i) -^|a)) is pure for any pure state |a)) = \a) (g) \a)* . 

(h) X'^X = Ijf»j2 = /„2 

(in) rank<I>(X) = 1 

Proof, (i) (ii). Condition (ii) immediately im- 
plies that |a')) = X\a)) is pure, since p[a'] — {{a'\a')) = 
((a|X^X|a)) = ~ 1- Conversely, as X can be 

represented by X = J2iLi Xi (81 X* , the quantum state 
I a')) = X\a)) always satisfies the following relation: 

((ay)) = ((a\X^X\a)) 

= (a\(ar Y^ixj <E> Xj){X, ® X*)\a)\a)* 

= ^|(a|4x,|a)p 

< J2{a\xlxy)(a\Xjx,\a) = 1. (16) 

Therefore, the condition ((a'\a')) = 1 imposes the equality 
relation in Eq. (|16p . Then, Xi\a) is parallel to Xj\a) for 
all and \a), indicating that Xi is independent of 

i. Thus, X takes the form X = X X*, where X is 
defined hy X :— ^/MXi. Consequently, using the trace- 
preserving condition X'^ X = lu, we arrive at X^ X = 

(ii) <^ (iii). First, we assume (iii). Then, ^{X) is 
written as <i>(X) = |a;))((a;| using a vector |a;)) & 1C(^H. 
Furthermore, as \x)) can be represented by |a;)) = {X (g) 
Ih)\Ih)) with a matrix X e £{n (g) /C), we have $(X) = 
{X(E) In)\In)) ((In\{X<» In)^, andthus X ^ X(g)X* from 
the definition of $. This directly yields X''X = /„2 due 
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to X'^ X = I-j-c- We next turn to the proof of (ii) (iii). 
Multiplying a pure state |a)) = |a) (g) |a)* e on both 
sides of /„2 = X^X where X — J2i ® obtain 

1 = {{a\a)) = {{a\X^X\a)) = J^K^lxj X,\a)\' 



< 



J2{a\XlX,\a){a\X^^X,\ 



Hence, from the same reason as in the proof of (i) 
(ii), X must be of the form X = X <Si X*. This im- 
pHes $(X) = (X (g) Ih)\Ih))({Ih\{X (g IhV and thus 
rank$(X) = 1. ■ 

Corollary 3. Suppose X G A'(7i,/C) satisfies 
rank$(X) = 1. Then, the nonzero eigenvalue of $(X) 
is given by n = dimTi. 

Proof. From the proof of Lemma 2, we have 

Tr$(X) = {{InliX^ X ^ In)\lH)) - {{Ih\Ih)) ^ n. U 

According to Lemma 2, the totality of quantum chan- 
nels that transform pure states in H to pure in /C is com- 
pletely characterized by the following non-convex set: 

Xi{n,ic) = e /:(7^®^/c^2) | i.ank$(x) = i, 

$(X) >0, {{I,c\X ^ {{In\ }. 



III. OPTIMAL ENCODER DESIGN AS A 
MATRIX OPTIMIZATION PROBLEM 

This section is devoted to rewrite the problem ([2]) as an 
encoder-optimization problem, which is further described 
as a matrix optimization problem using the notations 
introduced in Section II. 

First, let us fix the dimension of the codespace C to 
dimC = r and represent an element of C by \4>c) = E\(j)) 
with the input pure state \(j)) € C which contains all 
information of the sender. Here, E is the Kraus operator 
corresponding to the following encoding channel: 



E\(h) eCan. 



(17) 



We set H C"; then, i? is an n x r complex matrix 
satisfying E^E = Ir- In the example ([3]), |0) is a qubit 
and E = |00)(0| + |11)(1|, i.e., r = 2 and n 4. In 
terms of the above notations, the codespace-optimization 
problem ^ is written as 

P(yl)=max min P{A,EA(I))), 



P{A8M)) 



T,[A8{]<t.){<t>\f] 
Tv[A{E\c^)m^f] 



(18) 



Next, let us represent the problem using the matrix 
variable introduced in Section II. Since the encoding 
channel £ obviously preserves pure states, its matrix rep- 
resentation E is an element of A'i(C'',C"). Also, from 



Eq. mi), the input |0) takes the form \4>)) = in 
the extended space (C)**^. Hence, the output state of 
the encoder-error process is given by = 
where A E A'(C",C") is the matrix representation of 
the error channel A. Then, due to Eq. p5)) . the output 
purity is 

PiA,£, m = ((p'ip')> = {m\*E^A^AE\mr- 

Consequently, the max-min problem (jlSp is written as 



P(A) = max min {<j>\{<j>\* A^ AE\<I>)\<I>)* , (19) 

which is identical to the following "error-minimization" 
problem: 

min e, 

s.t. {(j)\{(j)\*E^A^AE\(j))\(j))* > 1 - e, V|0) G C\ 
E e A'l(C^C"), 

< e < 1. (20) 

Note that the optimal purity is related to the minimum 
error, Copt, by 

P{A) = 1 - Copt. 



IV. EXACT OPTIMAL SOLUTION TO THE 
PURITY-OPTIMIZATION PROBLEM 

In this section, we provide a systematic and powerful 
computational algorithm that exactly solves the purity- 
optimization problem when dimC = r = 2. The pro- 
posed algorithm can easily be extended to cover the gen- 
eral class of codespaces of dimension r > 3, in which case 
the suboptimized output purity gives a lower bound of 
the optimal purity P{A). This result will be discussed in 
Section V. 

The procedure to derive the algorithm consists of two 
stages. In the first stage, it will be proved that the first 
constraint in the problem (|^D)) : 



{cb\{cb\* E^ A^ AE\cbM* > 1 - e, V|(^) G 



(21) 



can be equivalently transformed to an LMI condition 
with respect to E, e, and an additional variable. In 
the second stage, we will consider a tractable rank- 
minimization problem of the matrix variable that is 
closely related to the original error-minimization prob- 
lem (PH)) . It will be then shown that, under a certain 
condition, the optimal solution of the rank- minimization 
problem coincides with that of the problem (|20p . 



A. The first stage: transformation of tlie constraint 

To simplify the exposition, we here assume that the in- 
put state |(/)) is a real-valued qubit, i.e., \(j>) = [xi,X2]^ € 
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(xf + xl = 1). The general qubit case |(/>) G will expressed as 



be discussed in Section IV-C using essentially the same 
idea presented here. 

Before considering the transformation of the constraint 
pip , let us further express it only in terms of real ma- 
trices. To this end, we define the following real matrix 
variable with the size 2n^ x 4: 



(22) 



Then, the output purity is expressed as 

p{A,£,m = {c^ME^pm\<j>), 

where P is a real positive semidefinite matrix defined by 

[ ^(A^A) n{A^A) 
Furthermore, we introduce a matrix (3^ 



B := 



1 

l/\/2 

l/\/2 

1 



(23) 



and define a vector 



\x))b :=BT|<^)|<^)=BT 



Xl 




Xl 


X2 




X2 



V^Xl 



X2 



Note that \x))b is normalized: b{{x\x))b = 1- As a result, 
from the relation 10) 10) — B\x))b, the constraint ([21]) is 



pix) B' E' PEB + {e ^ l)h\\x))B > 0, 

Vxi,X2eM. (24) 



We are now in the position to describe the transfor- 
mation. The constraint (|24p indicates that p{x) must 
be a real fourth-order nonnegative polynomial function 
with respect to the variables {xi,X2). This type of con- 
straint, i.e., the nonnegativity of a polynomial function, 
frequently appears in a wide variety of engineering prob- 
lems. In particular, the following Sum-of-Squares (SOS) 
characterization of non-negative polynomials, first stud- 
ied by David Hilbert more than a century ago, is a funda- 
mental question: When does a nonnegative polynomial 
p{x) have an SOS decomposition p{x) = J^i^ii^) ^^r 
some polynomials hi{x)l One of the well-known answers 
to the above question leads us to conclude that the non- 
negative polynomial p{x) must have an SOS decompo- 
sition, thereby Eq. (|24p is equivalently replaced by the 
following matrix inequality: 



B^E^PEB + (e - l)/3 + tS 

+ tIe^pet2 + T]E^PETi 

- TjE^PETi - tJe^PETs > 0, (25) 
where t g R is an additional optimization variable. The 
proof of Eq. (23) and the matrices Ti,T2,T3,Ti, and 
S are given in Appendix A. The inequality pS]) is trans- 
formed to 



TS+{e- l)/3 



" EB 


T 


" fc/2„2 - P 






ETi 






-P 




ET2 




-P 


fc/2„2 




ET3 






kl27i2 


p 


.ET4. 






p 


fc/2„2 





EB ' 




ETi 




ET2 




ET3 







k Ib'^e^eb + tIe^eTi 



tIe^et2 



tIe^et^ 



tIe^ETa 



>o, 



(26) 



where the blank spaces in the large matrix denote zero 
entries. The fixed scalar number fc > is selected such 
that 



fc/2„2 - P > 

is satisfied. This is equivalent to 



fc/2„2 P 
P fc^2n2 



> 0, 



fc/2„2 -P 

-P A:/2„2 



> 0. 



(27) 



(28) 



Then, due to the conditions ([27]) and ([SS]) . the large ma- 
trix in Eq. (|26p is positive definite. Moreover, we now 
see from Lemma 2 that the non-convex rank condition 
rank<i>(£^) = 1 is equivalent to E'^ E — I^, which leads to 

E'^E ^ ^{EY^{E) -f %Ey%E) = h. 

Thus, the last term in Eq. (f26|) is calculated to 

B^B + T[Ti + tJTs -I- tJTs + TlTi = 2/3. 
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used to transform Eq. p6l) to 
I 



p 



p 



EB 

ETi 
ET2 
ET3 
ETi 

(2fc + e-l)/3 



>0, 



(29) 



Finally, the Schur complement (see Appendix B) is 



(fcJ2„2 - P)-l 

-P fc/2„2 



B^E^ [tIe^ tIe^] 



which is obviously an LMI with respect to the variables 
E, e, and t. As a result, the original problem is equiva- 
lently written by 

min e, 

E,e,r 

s.t. (E,e,r)e7Vi, (30) 

where A/i is the following non-convex set: 

Ml := { (£;,£, r) I $(£;) > 0, {{I,,\E = {{hi 

LMI (IMl), < e < 1, rank$(£;) = 1 }. (31) 

B. The second stage: rank-minimization 

Let us consider a closely related problem 
min rank$(£J) + 7e, 

E,e,r 

S.t. (E,e,r)GA/', (32) 

where A/" is a convex set given by 

N:^ { {E,e,T) I $(£;) >0, {{I^\E ^ {{h\, 

LMI dSSl), < e < 1 }. (33) 

The tuning parameter 7 > gives the relative weight be- 
tween the two objectives rank<I>(£J) and e. This change of 
the problem is motivated by the fact that we can now ap- 
ply some known heuristic methods for rank minimization 
problems, one of which is discussed below. 

The minimization of the rank of a matrix subject to 
convex constraints is a ubiquitous problem in diverse ar- 
eas of engineering such as control theory, system identi- 
fication, statistics, signal processing, and computational 
geometry [2^. The general rank- minimization problem 

min rankA s.t. X e M and A > 0, 

where A > is the optimization matrix variable and A4 
is a convex set denoting the constraints, is computation- 
ally NP-hard, thus we need to rely on heuristics. The 
log-det heuristic introduced and discussed in [2^ [2^, [23| 



provides an attractive approach. The heuristic is de- 
scribed as follows: The function logdet(A + SI) is used 
as a smooth surrogate for rankA to yield 

min logdet(A -I- 51) s.t. X e M and A > 0, 

where (5 > is a small regularization constant, and 
can be chosen to be on the order of the eigenvalues we 
can consider as zero. Although the surrogate function 
logdet{X + SI) is not convex, it is smooth on the positive 
definite cone and can be minimized locally using any local 
minimization method; we here use iterative linearization. 
Let Aj denote the i-th iterate of the optimization vari- 
able A. The linearization of logdet(A + SI) around Xi 
is given by 

logdet(A + SI) = logdet(Ai + SI) 

+ Tr[(A, + «)-i(A-A,)],(34) 

where we have used the fact that VlogdetA = A^^ 
when A > 0. Hence, we can minimize logdet(A + SI) 
over the constraint set M by iteratively minimizing the 
local linearization ((34|) . This leads to 

X,+i = ar^minTr [(A, + SI)-^X] . 

The new optimal point is A^+i. Since the log-det func- 
tion is concave in A, at each iteration its value decreases, 
and the sequence of the function values generated con- 
verges to a local minimum of log det{X +SI) . This implies 
that the global optimal solution Aopt can be obtained by 
appropriately choosing an initial point Xq (see Fig.[T]). 

The above procedure is directly applicable to the case 
where the objective function is replaced by rankA -I- 
with e € [0, 1] an additional variable and 7 > a con- 
stant. Therefore, the rank-minimization problem (|32p is 
replaced by 

min log det ($(£;) + Shn) + 7e, 

E,e,T 

S.t. (E,e,r)eA/'. (35) 

The local or global optimal solution of this problem is 
obtained by solving the following iterative SDP: 

= argmin |Tr \{<^{Ei) + Shn)-^^{E)] + 75}. (36) 

(E,e,T)G7^/ I. -'J 
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Xf 



Xi' 



local minimum 
Xopt ■ global minimum 

FIG. 1: The log-det function and a convergence of the itera- 
tion variable Xi. 



Note that the convergence point of this algorithm is very 
sensitive to an initial point Eq. (We do not need to 
specify Eq ^nd tq , since and are not used to calculate 
{Ei+i,ei+i,Ti+i).) 

We now provide an important theorem that connects 
the replaced problem ((35| with the original problem ([30| . 

Theorem 4. If the local (global) optimal so- 
lution (-Eopt, Copt, Topt) of the problem ((35|) satisfies 
rank<l'(£Jopt) — 1, then it coincides with the local (global) 
optimal solution of the problem (jSO]) . 

Proof. The global optimal solution of ([55)1 satisfies 



logdet($(£;opt) +(5/2 n) ~r 7^opt 

< \ogdet{^{E)+Sl2n)+ie, 



(37) 



for all {E,e,T) e M. From the assumption and Corol- 
lary 3, the nonzero eigenvalue of <I>(£^opt) is dimC = 2, 
which yields det($(£;opt) + Shn) = (2 + S)S^''-\ Also, 
any E £ A'i(C2,C") satisfies det($(£;) + Shn) = [2 + 
(5)(5^"~^ As a resuh, Eq. reduces to 

Copt < e, V(£;,e,T) e TVi, 

where Mi is defined in Eq. pT|) . This implies that 
(£^opt, Copt, Topt) is indeed the global optimal solution of 
([30|) . We can prove the same fact for any local optimal 
solution by considering local regions of M and Mi . ■ 

Clearly, the above theorem can be extended to the gen- 
eral case of r. Therefore, the optimal (suboptimal) solu- 
tion of the original problem (^0)) is obtained by equiva- 
lently transforming (relaxing) it to a problem of the form 
(I30p and solving a related rank-minimization problem via 
the same heuristic. 

Finally, let us discuss choosing an initial point Eq 
of the algorithm ((36| such that the iteration vari- 
able (Ei,ei,Ti) converges to a local optimal solution 
(£^ AT, EAT, Tat) with rank$(£JAr) = 1. We here make the 
following observation; an initial point Eq that also sat- 
isfies rank<I>(£Jo) = 1 might be a good candidate for the 
above requirement to be satisfied. From the proof of 
Lemma 2, this implies 



Eq = Eq®E*, ElEo = h. 



(38) 



Actually, in many practical cases, we observe that an 
initial point of the form converges to a feasible local 
optimal solution. This fact will be seen in Section VI. 



C. Exact optimal encoder for general qubit inputs 

We here consider the general qubit input |0) = 
[e'" cos /3, sin /3] G (a, /? S M) and outline the equiva- 
lent transformation of the constraint (HJ). We first note 
that the input vector \4>)\4>)* is represented in terms of a 
monomial vector as follows: 



cos^ (3 
e'" sin fj cos (5 
e~'" sin /3 cos (3 

sin^/? 



U 



Xi + X2 

\f2xixz 

\/2x2X3 



= U\x))u, 



where the real variables xi,X2,X3 € R are defined as 

Xi = cos l3 cos a, X2 = cos /? sin a, x^ — sin /3, 
and [/ is a unitary matrix given by 
1 



U := 



I/V2 i/%/2 
1/V2 -i/x/2 



1 



Then, similar to the previous case, defining a 2n x 4 real 
matrix 



e' 



^{EU) 
^{EU) 



the output purity is expressed by 

P{A,£,m ^ u{{x\E^PE\x))u. 

Consequently, the original max-min problem is equal to 
the minimization of e S [0, 1] subject to the conditions 
E e A'i(C2,C") and 

p'ix) := u{{x\ [e"'PE + (e - 1)14] \x))u > 0, 

Since p'{x) is a fourth order homogeneous polynomial 
with respect to the three variables (xi, X2, 2:3), the 
Hubert's lemma (iii) in Eq. (jAl[) can be applied; the 
nonnegativity of p'{x) is equivalent to the condition 

p'{x) is an SOS with respect to (xi, 22, 2:3). 

Then, as the SOS decomposition ofp'{x) implies the ex- 
istence of a positive semidefinite matrix Q > satisfying 
p'{x) = u{{x\Q'\x))u, the matrix E^PE + (e - l)/4 is 
related to Q' by 



Q' = Y,T'J [e^PE + (e - 1)L 

i 

+ Y,T'r\E^PE+{e-l)l4 



rpff 

i 

T' 
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with certain matrices T'^,T^ , S'^, and additional scalar 
variables r^' G M. The above nonlinear matrix inequal- 
ity with respect to the variables E, e, and t- is further 
transformed to an LMI using the same technique shown 
in Section IV-A. As before, we then consider the problem 
of minimizing rank<i>(£J)+7e subject to the LMI obtained 
above and the linear constraints > 0, {{In\E = 

((/2I, and < e < 1. If the optimal solution of this prob- 
lem satisfies rank<l'(£^opt) = I7 then it is also the optimal 
solution of the original problem ([^01) with r = 2. 



V. SUBOPTIMAL SOLUTION IN HIGHER 
DIMENSIONAL CODE SPACE 

In the general case r > 3, nonnegativity of a homoge- 
neous polynomial no longer implies the existence of its 
SOS decomposition (this remarkable equivalence holds 
only in the cases (jAip ). However, the SOS characteriza- 
tion can still be used as a sufficient condition; that is, the 
first constraint in Eq. PP]) is relaxed to 

p"(x) := {{x\ [E^PE + (e - 1)/,.] \x)) 

is an SOS with respect to (xi, . . . ^X2r~i), (39) 



where E £ 



p2n xr 



is an appropriately defined real ma- 
trix variable that is linear to E, and \x)) € MJ" is an ap- 
propriately defined real monomial vector of xi, . . . , X2r-i- 
The SOS condition ((39l) equivalently leads to an LMI 
as seen before, and consequently, we have a problem of 
the form ([50)1 that can be tackled via the log-det heuris- 
tic. Note again that Eq. ([5^1) is only a sufficient condi- 
tion for the inequality p"{x) > to be satisfied for all 

(xi, . . . , X2r-i)- Thus, any feasible solution [E , e) satis- 
fying Eq. is included in the original set of solutions. 
Therefore, the suboptimal error computed from the re- 
laxed problem, Egub, is always bigger than or equal to the 
exact optimal error Copt- This indicates that the subop- 
timal output purity, Psub(-4) = 1 — Egub, gives a lower 
bound of the optimal purity: 



P{A) - 1 - eopt > 1 - es 



PsuhiA). 



An important fact to be noticed is that, as pointed out in 
[131 , the gap between the set of nonnegative polynomials 
and the set of polynomials with an SOS decomposition is 
considered to be small in a practical situation. Hence, we 
expect that Psuh{A) is a good approximation to P{A). 



VI. EXAMPLES 

A. The bit-flip channel 

The quantum bit-flip channel with flipping probability 
p is given by 

5(C^) 3 Tip = pa^pTx +qp £ 5(C^), 



where p+q = 1 and (Tx = |0)(1| + |1)(0|. We here consider 
the double bit-flip channel Am = Tf^: 

4 

S{C^) 3p^p'^ Amp = ^-P^\ € 

1=1 

Ai=pax® <7x, A2 = ^/pqcTx <E) I2, 
A3 = h^fJx, Ai = qh® h- 



The matrix form of the double bit-flip channel. Am 
J2i Ai ® A* e A'(C'', C), is represented by 





- qAi 


^/pqAs 


^/pqA2 


pAi 


Am = 




qAi 


pAi 


^fvqA2 


y/pqA2 


pAi 


qAi 






. pAi 


y/pqA2 




qA^ 



In particular, we set p ~ 0.1; then, for example, k — 2 
satisfies the condition (f27|) : kl^2 — P > 0- 

We here assume that the input is a real- valued qubit: 
|(/)) e R'^. Then, an exact local or global optimal en- 
coder E„ 



■'opt 



e ^"1(1 



) is computed by the algorithm 



([36|) under the condition rank$(£Jopt) = 1- A strong 
convergence property of Ei is observed when the SDP 
parameters are set to S = 0.01 and 7 = 15. We usually 
need 90 iterations of the SDP; hence we denote the con- 
vergence point by (iJgo, £90, T90). Note again that Ego 



Eg 



Eon due to the rank 



must be of the form Ego 
condition rank$(£J9o) = 1. Regarding the initial point 
£^0? we follow the idea mentioned in the last paragraph 
of Section IV-B and examine some initial points of the 
form ((38)) to find the global optimal solution. 

First, we randomly choose two initial points as E^ = 



Eq^ (g) Eq'* {j = 1,2), where the Kraus operators Eq 
and eI^^ are given by 



(1) 



E, 



(1) 



2 
V2 -V6 
V3 1 

1 V3 



e: 



(2) 



^/2 

V3 -76 

1 \/2 

2 ^/2 



respectively. Then, the corresponding convergence points 

(i = l,2). 



(j) 



are respectively given by Ego 
where 



E, 



U) 

90 



E, 



U)* 



i ^90 



77.(2) 1 

0.5308 -0.4672 

0.5308 -0.4672 

0.4672 0.5308 

0.4672 0.5308 



0.5274 -0.4710 

0.5274 -0.4710 

0.4710 0.5274 

0.4710 0.5274 



In both cases, the convergence value of the error is given 



by £90 = 0.18. In view of the structure of e'^q and K 
we expect that the encoder E^^^: 



_g.(a) ^ £;(a)^£;(a)* £;(a) ^ 



1 

3 



90 

cos a — sin a 
cos a — sin a 
sin a cos a 
sin a cos a 



,-.(2) 
^90 ' 



(40) 
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will be a local optimal solution and provide a local mini- 
mum of the error, e*^^) — 0.18, for all a € [0, 27r). In fact, 
for the input \(p) = [xi, X2]^ — [cos 95, sin 1^9]''' and the 
encoder £'-^\ the output purity is reduced to 

2a)]', 



P(Af,f<"\|0)) = 1-2to[cos(2(p 
which takes the minimum value 



P. 



(a) 



minP(Af,f("M<^» 



l-2pq = 0.82. 



Hence, as expected above, the local minimum of the error 
is e('^) = 1 - 0.82 = 0.18. This result clarifies that the 
optimal encoder depends on the worst-case input as a = 
— "y^worst + nT:/2, where n is any integer. As a summary, 
the encoder 



£(a) 



1 
72 



xi cos a — X2 sma 
xi cos a — X2 sin a 
xi sin a + 2:2 cos a 
xi sin a -I- 2:2 cos a 



is locally optimal for all a G [0, 2tt). 



We next try following two initial points, jj/q 



E 



U), 



E^^>* {j = 3,4), where 



p(3) _ 
Eq - 



2 
V2 V6 
V3 -1 
-1 V3 



^0 ^ 



a/3 
V2 

2 

1 



-1 


-a/3 



Then, the corresponding convergence points are respec- 
tively given by 



90 



E, 



U) 

90 



4o* (j = 3,4) where 



r p(3) p(4) . 

1 ^90 J ^90 / 

0.6935 
0.1377 
0.6935 
0.1377 



-0.1377 
0.6935 

-0.1377 
0.6935 



0.4636 
0.5341 
0.5341 
0.4636 



-0.5341 
0.4636 
0.4636 

-0.5341 



Although they have a similar structure, there is a large 
gap between the corresponding convergence values of the 
error e: 



,(3) _ 



0.18, 



^90 



= 0.2952. 



The structure of E^q and Eqq suggests that the encoders 
Eit^) = ^(m) (g, ^(m)* = b, c) with 

{ e'^^\ S^'^) } 
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V2 



cos a — sm a 

sin a cos a 

cos a — sin a 

sin a cos a 



1 
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cos a — sm a 

sin a cos a 

sin a cos a 

cos a — sin a 



are locally optimal for all a G [0, 27r). Actually, the out- 
put purity with the above encoders and the input 
10) = [cosifi, sin ip]^ are respectively calculated as 



P(Af,f('^\|0)) = i 
P(At,f(^\|0)) = i- 



-2pq[cos(2(^-|-2a)]2, 
4pq{p^ + q^)[cos{2Lp + 2a)]^ . 



Thus, their minimum values are 

P!nl = ^-2pq = 0.82, 

+ q^) = 0.7048, 



P^") = 1 

min 



irrespective of a. The minimums are attained when 
cos(2iy9 + 2a) — ±1, as in the case of S'-^K We also see 
the following inequality: 



p(b) _ p(c) 
min min 



2^(1 - 2pf > 0. 



Therefore, the encoders S'^'^^ and S^^'' achieve the same 
local minimum of the error, whereas S^'^'^ is inferior to 
those channels for all p. 

Combining the entire set of investigations presented 
above with other numerical results that were omitted for 
brevity, we maintain that Copt = 0.18 is the global mini- 
mum and that the optimal purity is thus given by 



P{A 



hf) 



1 - eopt = 0.82. 



The solutions S'-'''^ and S'^^'^ are typical optimal encoders 
that yield the above optimal purity. 

Remark 1. In the Kraus representation, the 
output state is given by p' = J2i^iE\(j)){(j)\E^Al. 
Intuitively, in order for the purity of p' to have a large 
value, the encoder E should be chosen so that the 
vectors {AiE\(f>)} are close to each other. Actually, if all 
of them are parallel, the output state is pure. In this 
sense, f^'*' is a physically reasonable encoder because 
the vectors AiE^^'^cj)) and As^J^'^^l^) are parallel to 
A2E'''^^ and A4E'-''^\4i) , respectively. The encoders 
E^^') and i?*^'^-' also satisfy such relations. In contrast, 
if we choose E = |00)(0| + |11)(1| in Eq. the four 
vectors AiE\(l)) (i = 1, . . . ,4) differ from each other and 
span a linear space of dimension 4. This is indeed a bad 
encoder since the minimum output purity in this case is 
calculated as p[p'] — {p^ + q^Y ~ 0.67, which is clearly 
less than the optimal purity P(Af) = 0.82. 

Remark 2. We again maintain that Eq satisfying 
rank<I>(£Jo) = 1 is a good initial point. Actually, within 
our investigation, we have observed that such an initial 
point always converges to a rank-one solution by appro- 
priately choosing the SDP parameters 5 and 7. However, 
for initial points with the rank more than one, it is easy 
to find a bad example of £^0 such that rank$(£Jgo) — 1 
is not achieved for any 5 and 7. For instance, if we 
choose <^{Eq) = (l/4)/8, then Ei always converges to 
a solution satisfying rank<I>(£J9o) = 2. Another rea- 
son of the emphasis is based on the following observa- 
tion. Once we obtain a rank-one solution using an ini- 
tial point with the rank more than one, then we always 
have found a rank-one initial point that converges to the 
same solution. In other words, it is considered that any 
rank-one solution is available by choosing a rank-one ini- 
tial point appropriately. For example, Ei starting from 



Eo = 0.5£;, 



(1) 



o.3£;: 



(2) 



one solution of the form 



(3) 

Q.2Eq converges into a rank- 
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B. The amplitude damping channel 

The amplitude damping channel describes the dissipa- 
tion of a quantum state into equilibrium due to coupling 
with its environment. The Kraus representation of the 
channel is given by 

'■)^ P^r2p^ HipHl + H2pHl e S{C^), 



where 



Hi = 



1 

VP 



H2 = 









The parameter p e (0, 1) represents the rate of dissipa- 
tion. We consider the double amplitude damping channel 

4 

1=1 

Ai iJi (g) iJi, A2 = i/i (g) H2, 
A:i^H2® Hi, Ai^H2® H2. 

The matrix form of the channel, Aad = X^i Ai ® A* e 
A'(C4,C4), is given by 



A 



ad 



Ai 


VI - PA2 


%/l - pAs 


{l-p)Ai 1 


O4 


VPAi 


O4 


Vp(1-p)A3 


O4 


O4 




Vp(1-p)A2 


O4 


O4 


O4 


Ml 



where O4 denotes the 4x4 zero matrix. In particular, we 
consider the case of p = 0.1 and set fc = 4, which leads 
to fc/32 - P > 0. 

Our goal is to obtain the optimal encoder under the 
condition 10) e in which case Eopt e Xi{R^,C'^). 
The iteration variable Ei of the algorithm is initial- 
ized to Eq of the form ((38)) . and the SDP parameters 
are set to 6 — 0.01 and 7 = 6.1. In order to find a 
rank-one convergence point, we usually need 500 itera- 
tions of the SDP; we thus denote the convergence point 
by (-E500, £500, Tsoo)- 

First, let us take the initial points i^Q^-* and -Eq^\ which 
have appeared in the bit-flip channel case. Then, the cor- 
responding convergence points are respectively given by 



E 



U) 

500 



-^500 ® -^500* (i — 5j 6) with the Kraus operators 



r p(5) I 

t -^500' ^500 J 



(6) 



0.6555 0.2503 

0.4174 -0.9045 

0.5741 0.2822 

0.2583 0.1989 



0.5803 -0.2273 

0.4668 0.8745 

0.5568 -0.2843 

0.3679 -0.3209 



In both cases, the convergence value of the error is given 
by £500 — 0.18. Unlike the case of bit-flip channel, the 
above solutions do not have a simple structure of the ma- 
trix entries, which is highly important for a physical real- 
ization of encoding process. To obtain a simple solution. 



let us carry out the algorithm with an initial point that 
has a specific matrix form itself. As a typical example, 
we consider the following initial point: 



.(d) 



E, 



(d) 



77.(d)* p(d) 
^0 ' ^0 



cos a 

cos/3 

sin a 

sin/3 



Then, for any a £ [0, 27r) and (3 S (0, n/2), Ei converges 
to 



ir(d) _ p(d) „ p(d)* rpy 
^500 ~ ^500 ^ ^500 ' -^500 



r(d) 



COS a 

1 

sin a 





with egQQ — 0.18. This encoder is locally optimal for 
all a £ [0, 27r). Actually, the output purity (fT8)) for 
the encoder-error process And^^oo with the input 10) = 
[xi, X2]^ is calculated as 

:-(d) _ 1 r,„/l „\/^2„;„2, 



P(Ad,f^oo, m = 1 - 2pil-p){xi sin^ 



and thus, its minimum is P 



(d) 



1 - 2p(l - p) 



(41) 
0.82 



when 10) = [0, 1]^ irrespective of a. 

We also observe the following similar convergence: 

cos a ■ 

sin a 

1 





^0 - 



cos a 

sin a 

cos/3 

sin/3 



^500 ^ 



with e^ol, 0.18 for all a £ [0, 27r) and /3 £ (0, n/2). The 
output purity P^AndT^^oo, 1*^)) same form as Eq. 

(PTjl . thus the encoder ^^qq is also locally optimal for all 
a £ [0,27r). 

Finally, let us choose an initial point of the form 



'E, 



(f). 



^0 — 



cos a 

cos (3 

sin/3 

sin a 



(42) 



We then observe a somewhat complicated convergence 
depending on {a, (3) as follows. When a takes a small 
number, e.g., a — 0.2 (any /3 can be taken), the algorithm 
does not cause a variation in Ei, and only changes into 
0.18. That is, we obtain the local optimal solution 



E 



500 



e: 



(fi) 

500 



'e: 



(fi). 

500 



E. 



(fi) 
■500 



cos a 

cos/3 

sin/3 

sin a 



On the other hand, when a ~ 7r/2, another type of con- 
vergence occurs. For example when choosing a — 1.3, Ei 
converges to 



E 



ih) _ p(f2) 



500 



E, 



(f2). 



e: 



(f2) _ 

500 " 



0.6893 

cos/3 

sin /3 
0.7245 
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with e 



(f2) 

500 



0.18. To further understand this complex 



structure of the solution, we provide an analytical in- 
vestigation of the output purity P{Aiid,£o \ \4>)) in Ap- 
pendix C. However, we reemphasize that a lucid advan- 
tage of our method to search an optimal solution is that 
it does not require any analytic examination on the max- 
min optimization problem of the output purity, which is 
in general extremely hard. 

Based on the above investigations, we maintain that 
£opt = 0-18 is the global minimum and that S^^qq (/i = 
d, e, fi,f2) are the optimal encoders. Therefore, the opti- 
mal purity is given by 



■ad) 



1 



topt 



0.82. 



VII. CONCLUSION 

In this paper, we presented a tractable computational 
algorithm for designing a quantum encoder that maxi- 
mizes the worst-case output purity of a given decohering 
channel over all possible pure inputs. We cast the prob- 
lem as a max-min optimization problem (minimization 
over all pure inputs, and maximization over all pure state 
preserving encoders) . Although this problem is computa- 
tionally very hard to solve due to the non-convexity prop- 
erty, our algorithm computes the exact optimal solution 
for codespace of dimension two. Moreover, we showed an 
extended version of the above algorithm that computes 
a lower bound of the optimal purity for the general class 
of codespaces. 

We believe that the proposed computational approach 
provides a powerful method that is also applicable to 
other problems in quantum encoding and fault-tolerant 
quantum information transmissions. For example, fol- 
lowing the same techniques presented in this paper, we 
can prove that a quantum error correction problem with 
the minimum fidelity criterion considered in [Tsl is 
transformed or relaxed to a convex optimization prob- 
lem systematically; we are then able to obtain the opti- 
mal or suboptimal solution using SDP. This result will 
be reported soon. 
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APPENDIX A: PROOF OF EQ. ((25l) 

Let us consider a real polynomial function p{x) in n 
variables x = (xi, . . . , Xn) of the form: 



Cfe e 



where the sum is over n-tuples k = (fci, . . . , A:„) satisfying 
~ This function is called the homogeneous 
polynomial of degree m in n variables. A homogeneous 
polynomial satisfies p{Xxi, . . . , \xn) = A™p(xi, . . . , Xn). 
We now state the famous Hilbert's theorem. Let Pn,m 
be the set of nonnegative homogeneous polynomials of 
degree m in n variables. Let S„^m be the set of homoge- 
neous polynomials p{x) that has an SOS decomposition 
p{x) = hi{x)'^ , where hi{x) are homogeneous polyno- 
mials of degree m/2. Then, Pn,m — '^n,m holds only in 
the following cases: 

(i) n = 2 (ii) m = 2 (iii) n = 3, to = 4. (Al) 

For more detailed description on this problem, see (3l| . 
Now, Eq. (I24p has the following form: 



:){x) — [xi V2xiX2 x\]H 



^/2xl 



X2 



x\ 



> 0, yxi,x2 e 



(A2) 

where H = (hij) is a real 3x3 symmetric matrix. The 
function p{x) is a homogeneous polynomial with respect 
to two variables xi and X2 (and degree m = 4). There- 
fore, from the Hilbert's formula (i) in Eq. (|A1[) . the con- 
straint (|A2[) is equivalent to the condition 

p{x) is an SOS with respect to xi and X2- 

Moreover, it can be shown that the existence of an SOS 
decomposition is equivalent to the existence of a positive 
semidefinite matrix Q = (qij) > such that 



p{x) — z{x)'^ Qz{x), 



(A3) 



where z{x) is a vector of monomials of degree equal 
to deg(p)/2 = 2. Comparing Eq. with (M\i . 

we set z{x) = [x^, ^/2xiX2, a;!]""" 
z{x)'^ H z{x) = z{xyQz{x) yields 

hii ^ qii, hi2 = gi2, his + 

h23 — 923, ^33 = 933, 

which leads to 



. Then, the equality 

^22 = 9l3 + 922, 



Q = 

As a result, Eq. 
inequality: 



hii 
hi2 

913 



hi2 

- hi3 

h23 



913 



913 
h23 



> 0. 



2]) is equivalent to the following matrix 



/ill hi2 

hl2 h22 + hi3 /l23 
ft.23 ^33 





" 





1 " 


+ r 





-1 







1 









> 0, 
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where r := 513 G R is an additional optimization vari- 
able. The above inequality can be expressed as 

H + S'lHS2 + SlHSi - SjHSi - SJHSs + tS>0, 

where 



Hence, the following relation holds: 



Si = 



1/V2 






S2 = 





' 1 





0" 




" 





0' 


83 = 











, S'4 — 






























1 






1/V2 



s = 



From the above discussion, the constraint 
lently transformed to 



1 
0-10 

1 



IS equiva- 



B^E^PEB + (e - l)/3 + 


tS 




+ sl 


B^E^PEB + (e 




^2 


+ sl 


B^E^PEB + (e 


- l)/3 


^1 


-si 


B^E^PEB + (e 


-l)/3 


^4 


-si 


B^E^PEB + (e 


-l)/3 


^3>0 



As SIS2 = O and SIS4 = O, we obtain Eq. ({251) : 
B^E^PEB + (e - l)/3 + tS 

+ tJe^pet2 + tIe^peTi 

- TlE^PETi - tIe^PETs > 0, 
where the matrices Ti {i — 1, ... ,4) are defined as 
Ti :— BSi, T2 :— BS2, T3 :— BS3, T4 := BS^. 

APPENDIX B: THE SCHUR COMPLEMENT 

The Schur complement is a powerful tool that trans- 
forms a convex but nonlinear constraint with respect to 
matrix variables into an equivalent LMI. Its derivation is 
very easy; assuming ^ > 0, we have a matrix equation of 
the form: 

/ -A-^B 
O I 



I 


0' 




■ A B 








St c 




■ A 









C 


-B^A- 



A B 
fit c 



> ^ 



A>Q 

C-B^A-^B > 0. 



This is termed the Schur complement. In order to see 
the usefulness, let us consider a nonlinear constraint of a 
matrix variable X: I — X^'X > 0. The Schur complement 
states that the constraint is equivalent to 



/ X 
At / 



>0, 



which is obviously an LMI. 

APPENDIX C: AN ANALYTIC INVESTIGATION 
OF THE PURITY-OPTIMIZATION PROBLEM 



We here give an observation on the purity-optimization 
problem where the error channel is And and the encoder is 
= e'^^pE'^P* with e'^^^ given in Eq. The out- 

put purity P{AnA,£^i\m - Tr[Adf,$'\|</')(0|)'] with 
the input \(t>) = [xi, X2Y '= is then calculated to 



P{a, /3, Xi) = 1 — 2pq (1 + sin 2q; sin 2/3 — 2pqs\n^ a)x\ 
— (1 -|- cos 2a + sin 2a sin 

First, let us consider the case where a takes a small num- 
ber. Especially when a = 0, P(0,/3, a;i) = \ — 2pq{x\~l)'^ 
is a concave function with respect to xi. Thus, the min- 
imum is given by Pmin — -P(0, /?, 0) = 1 — 2pq at xi — 0. 
This fact is still true for a « 0; the function P(a,/3,a;i) 
is concave and takes the minimum 1 — 2pq at xi = 
without respect to the values of a and (3. This is the rea- 
son why a and (3 do not have specific optimal values and 
the iterative SDP initialized with a « does not renew 
these parameters. On the other hand, when a = 7r/2, 
the output purity becomes 



P(7r/2, /3, xi) = 1 - 2pq (/ + q^)x\ + 1 



which obviously takes the minimum at xi = 1. Moreover, 
for a ~ tt/2 the function P{a, P, xi) is still concave and 
takes the minimum P{a, = 1 + Ap^q^ siv? Q;(sin^ a — 
1/pq). Unlike the case of a « 0, this function must be 
further maximized with respect to a. For this reason, 
there is a specific optimal value of a, whereas /3 does not 
affect the optimality. 
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